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CHARACTERISATION AND REPRESENTATION OF NON-DISSIPATIVE 
ELECTROMAGNETIC MEDIUM WITH A DOUBLE LIGHT CONE 

MATIAS F. DAHL 



^SJ ' Abstract. We study Maxwell's equations on a 4-manifold A^ with a medium that is non- 

dissipative and has a linear and pointwise response. In this setting, the medium can be 
^^ , represented by a suitable (j) -tensor on the 4-manifold A'^. Moreover, in each cotangent 

Cn ' space on N, the medium defines a Fresnel surface. Essentially, the Fresnel surface is a 

$_j I tensorial analogue of the dispersion equation that describes the response of the medium 

for signals in the geometric optics limit. For example, in isotropic medium the Fresnel 
surface is at each point a Lorentz light cone. In a recent paper, I. Lindell, A. Favaro and 
L. Bergamin introduced a condition that constrains the polarisation for plane waves. In this 
paper we show (under suitable assumptions) that a slight strengthening of this condition 
gives a pointwise characterisation of all medium tensors for which the Fresnel surface is 
the union of two distinct Lorentz null cones. This is for example the behaviour of uniaxial 

r^ ' medium like calcite. Moreover, using the representation formulas from Lindell et al. we 

CIh! obtain a closed form representation formula that pointwise parameterises all medium ten- 

sors for which the Fresnel surface is the union of two distinct Lorentz null cones. Both 
the characterisation and the representation formula are tensorial and do not depend on local 

Cy \ coordinates. 

1. Introduction 

vQ I We will study the pre-metric Maxwell's equations, where Maxwell's equations are written 

(^ • on a 4-manifold N and the electromagnetic medium is described by a suitable antisym- 

cn , metric (2) -tensor k on N that pointwise is determined by 36 real parameters. In each 

. I cotangent space on A^, the electromagnetic medium determines a fourth order polynomial 

^ ■ surface called the Fresnel surface that can be seen as a tensorial analogue of the dispersion 

b:^ ! equation. The Fresnel surface describes the response of the medium to signals in the geo- 

metric optics hmit [.OFROOl lRub02l IHO03[ lPSW09. .RRS 1 1 .1 . In this work we will assume 
that the medium is skewon-free. Then there are only 21 free parameters and such medium 
. , models non-dissipative medium. For example, under suitable assumptions the skewon-free 

r> ■ assumption will imply that Poynting's theorem holds I IHO031 iDahlO I. On an orientable 

C^ . manifold one can show that invertible skewon-free (2) -tensors are in one-to-one correspon- 

dence with area metric. By an area metric, we here mean a (^) -tensor on A^ that defines a 
symmetric non-degenerate inner product for bivectors. Area metrics appear when studying 
the propagation of a photon in a vacuum with a first order correction from quantum electro- 
dynamics [DH80, SWWIO]. The Einstein field equations have also been generalised into 
equations where the unknown field is an area metric IIPSW071 . For further examples, see 
I1PSW091ISWWI0I . 

We know that in isotropic medium like vacuum, the Fresnel surface is a Lorentz null 
cone at each point in N. That is, Lorentz geometry describes the propagation of light in 
isotropic medium. Conversely, it was conjectured in 1999 by Y. Obukhov and F. Hehl 
IIOH991IQFR00II that isotropic medium is the only (non-dissipative and axion-free) medium 
where the Fresnel surface is a Lorentz null cone. This was partially proven already in 
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liOFROOl . However, the full conjecture was only established in BFBllll by A. Favaro and 
L. Bergamin. For an alternative proof, see lIDahllal and for further discussions and related 
results, see IIQR021 lHO03l ILH041 Ilti05l and Section [U below. 

Since the Fresnel surface is a 4th order polynomial surface, the Fresnel surface can also de- 
compose into the union of two distinct Lorentz null cones. For example, this is the case in 
uniaxial medium like calcite (CaCOa) IBW99. Section 15.3]. In such medium, the propaga- 
tion properties of the medium does not only depend on direction, but also on the polarisation 
of the wave. In uniaxial medium, there are two eigenpolarisations and one null cone for each 
polarisation. In consequence, there is one Fermat's principle for each polarisation [|PSW09J . 
This is the the source for the physical phenomenon of double refraction. 

We know that uniaxial medium is an example of medium with two distinct null cones. 
A natural next task is to understand the structure of all medium tensors with this prop- 
erty. This is the main result in [Dahllbl, which gives the complete local description of 
all non-dissipative medium tensors for which the Fresnel surface is a double light cone (up 
to suitable assumptions). The importance of this result is that it shows that are three and 
only three medium classes with this behaviour. Moreover, the theorem gives explicit coor- 
dinate expressions for each medium class. The first medium class is a slight generalisation 
of uniaxial medium. The second class seems to be a new class of mediums. The last class 
seems to be unphysical; heuristic arguments and preliminary numerical tests suggest that 
Maxwell's equations are not hyperbolic in that class HDahllbi In the below, this result is 
summarised in Theorem [ 



The main contribution of this paper is Theorem 15. II Under suitable assumptions, this theo- 
rem gives a tensorial characterisation (condition |f//j| in Theorem 15. II ) of all non-dissipative 
medium tensors for which the Fresnel surface is two distinct light cones. In a suitable limit, 
the condition also reduces to the closure condition ¥? = —Aid for a A > that charac- 
terises medium with a single light cone IIHQ031 . Moreover, in Theorem 15. II we give a ten- 
sorial representation formula (equation (l64l ) that parameterises all non-dissipative medium 
tensors with two distinct light cones. Both the characterisation and representation formula 
are pointwise results. 

The background and motivation for Theorem 15.11 comes from a recent paper by I. Lin- 
dell, A. Favaro and L. Bergamin [LBF12|. In Section |4] we will briefly summarise some 
of the results from IILBF12L In this paper, the authors introduces a second order polyno- 
mial condition on the medium tensor (equation (l54l ) in the below). Equation (l54l ) is derived 
from a constraint on polarisation of plane waves, and in I1LBF12I1 it is shown that whenever 
condition (l54l ) is satisfied (plus some additional assumptions), the Fresnel surface always 
factorises into two second order surfaces. In Section l43] we will further motivate that equa- 
tion (l54l ) is in fact a general factorisability condition for the Fresnel surface. At first this 
might seem unexpected since equation (l54l ) was initially derived from a constraint on polar- 
isation, yet it is able to constrain the behaviour of signal speed. However, the explanation is 
that for electromagnetic waves, polarisation and signal speed are not independent properties 
but tied together. In Theorem 15. II condition jf/ZJl is a slight strengthening of equation (l54l ). 
Also, representation formula (l64T l in Theorem 15. H is adapted from |LBF12| and constitute 
a subclass of generalised Q-medium introduced by I. Lindell and H. Wallen in IILW02II . A 
further technical discussion on Theorem 15. H is given in the end of Section [5] 

Some of the computations in the paper rely on computer algebra. For further information 
about the Mathematica notebooks for these computations, please see the author's homepage. 
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2. Preliminaries 

By a manifold N we mean a second countable topological Hausdorff space that is locally 
homeomorphic to M" with C°°-smooth transition maps. All objects are assumed to be 
smooth where defined. Let TN and T*N be the tangent and cotangent bundles, respec- 
tively. For A: > 1, let Q^{N) be antisymmetric tensor fields with k lower indices (that is, 
/c-forms). Similarly, let rjfc(A^) be antisymmetric tensor fields with k upper indices. More- 
over, let i7S(A^) = Q'^{N) (g) ^2{N). Let also C°°{N) be the set of scalar functions (that 
is, (q) -tensors). The Einstein summing convention is used throughout. When writing ten- 
sors in local coordinates we assume that the components satisfy the same symmetries as the 
tensor. 

2.1. Twisted tensors. If N is not orientable we will also need twisted tensors IIHO031 
Section A.2.6]. We will denoted these by a tilde over the tensor space. For example, by 
Q'^{N) we denote the space of twisted 2-forms. If G G Q'^{N) then in each coordinate 
chart ([/, X*), G is determined by a usual 2-form G\u G i7^(C/) and on overlapping charts 
([/, x*) and ([/, x*), forms G\u and G\fj satisfy the transformation rule 

(1) G\^ = sgndet ^-^JG\u, 

where sgn : R — )• M is the sign function, sgn x = x/\x\ ior x ^ and sgn x = for x = 0. 
If locally 

(2) G\u = -Gijdx^Adx^, G\fj = -Gijd¥ A dx\ 
then equation ([T]) implies that components Gij and Gij transform as 

When the chart is clear from context, we will simply write G = ^Gijdx^ A dx^ . Similarly, 
if K G Q'^2iN) then in each chart k is represented by a k|(/ G $7^2 (f^) and locally 

(4) K = -Kiidx'^Adx'^—^A^ 

for suitable components Krs- Moreover, if K^i and Kr4 are components for k in overlapping 
charts (C/, x*) and (t/, x*) then we obtain the transformation rule 

. _ . (dx-^\ p,dx^dx;^dx^9x^ 

Compositions involving twisted tensors are computed in the natural way by composing 
local tensors. For example, if k, ry G Q^2(-^) their composition defines an element ko r] £ 
Q'^2{N) and if k, t] and k o ry are written as in equation (HJ then 

(6) i^orjti = \trf^. 

If M is orientable, then twisted tensors coincide with their normal (or untwisted) coun- 
terparts. For example, if M is orientable, equation ([5]) implies that (7^2 (-^) = ^'^2{N). 
There are also other way to define twisted forms. Equation ([1} coincides with definition 
of a pseudo-form in IIFra041 . For a global definition of twisted forms using the orientation 
bundle, see lAMROli Supplement 7.2A]. 
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2.2. Tensor densities. In addition to tensors and twisted tensors, we will need tensor den- 
sities and twisted tensor densities. A (^\-tensor density of weight w ^ Zona, manifold A^ is 

determined by components T^ \"'f^ ^ in each chart ([/, x*), and on overlapping charts ([/, x*) 
and {UjX"^) we have the transformation rule | Spi99) , 



A twisted (^^^-tensor density of weight u; G Z on A^ is defined in the same way, but with an 
additional sgndet ( |^ j factor in the transformation rule as in equations © and (|5). 

The Levi-Civita permutation symbols are denoted by Siju and e*-'^'. Even if these coincide 
as combinatorial functions so that eijki = e^^^\ they are also different as they globally 
define different objects on a manifold. Namely, if eijuue"^^^^ and eijku ?"'^' are defined on 
overlapping coordinate charts ([/, x*) and ([7,5;*), respectively, then 

_ / 5x* \ dx"^ dx"^ dx^ dx^ 

(') £abcd = det ( -j^— T ) Epqrs 



dxi J ^'^^^ dx"- dx^ dx^ dx'^ ' 
^bcd _ .. (dxl_\ dx^&x^dx^dx^ 



\ dx^ J dxP dxi dx'' dx^ 

That is, Eijki defines a (^) -tensor density of weight — 1 on A^ and e*-''^' defines a (q) -tensor 
density of weight 1. For future reference, let us note that 

(9) ^^""'er-sii = 4<5f,5;] , e™^%ri,fc = 'i^S^^ , 

where (5* is the Kronecker delta symbol and brackets [ii . . . ip] indicate that indices ii, . . . ,ip 
are antisymmetrised with scaling I /pi. 

2.3. Maxwell's equations on a 4-manifold. On a 4-manifold A^, the premetric Maxwell's 
equations read 

(10) dF = 0, 

(11) dG = J, 

(12) G = k{F). 

where d is the exterior derivative, F G ^^(A^), G G fl'^{N), J G h'^{N) and k G Q'^2{N). 
Here, F, G, are called the electromagnetic field variables, J describes the electromag- 
netic sources, tensor k models the electromagnetic medium and equation (IT2l) is known 
as the constitutive equation. In local coordinates, equations (IT0l)- (IT2l) reduce to the usual 
Maxwell's equations. For a systematic treatment, see IIRub02 . HO03|. 

If locally F = ^Fijdx^ A dx^ , G = ^Gijdx'' A dx^ and k is written as in equation ^ then 
constitutive equation (fT2l ) is equivalent with 

(13) Gij = -KijFab- 



Thus equation jTzh models electromagnetic medium with a linear and pointwise response. 

Suppose K G Q'^2{N) and suppose {U,x^) is a chart. Then the local representation of k 
in equation ^ defines a pointwise linear map Q'^{U) — ;■ ^"^{U). In U we can therefore 
represent k by a smoothly varying 6x6 matrix. To do this, let O be the ordered set of index 
pairs {01, 02, 03, 23, 31, 12}, and if J G O, let dx-^ = dx'^^ A dx'^\ where Ji and J2 are 
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the individual indices for J. Say, if J = 31 then dx^^ = dx^ A dx^. Then a basis for il^(C/) 
is given by {dx'^ : J G O}, that is, 

(14) {dx^ A dx'^ , dx^ A dx^ , dx^ Adx^,dx^ Adx^,dx^ Adx^,dx^ Adx^}. 
This choice of basis follows IIHO03[ Section A. 1.10]. By equation dUl it follows that 

(15) K{dx-^) = ^KJdx^, JeO, 

leo 



where kJ = Kj^j^- Let b be the natural bijection b: O — )• {!,... ,6}. Then we identify 
coefficients {kJ : I , J £ 0} for k with the smoothly varying 6x6 matrix P = {nj)ij 
defined as nj = Pb{i)b{j) for I,J€0. 

Suppose P = {k'J)ij and P = (k/)/j are smoothly varying 6x6 matrices that represent 
tensor k in overlapping charts {U, x*) and ([/, 5;*). Then equation © is equivalent with 



_ 7 , / 5x* \ v^ dx^ T dx'^ ^ ^ ^ 



where 



dx'^ dx"^^ dx'^'^ dx-^'^ dx'^^ 



dx^ 9x^1 0x^2 Qxh Qxh ' 

and ^-[ is defined similarly by exchanging x and x. For matrices T = (^j)/j and 5 = 

{^t)ij, we have T = S~^, whence equation i^ is further equivalent with the matrix 
equation 

(17) P = sgndetf^j TPT-\ 

In a chart ([/, x*), we define trace k: [/ — )• M and det k: [/ — )• M as the trace and deter- 
minant of the pointwise linear map Q'^{U) — )• 0^(C/). When P is as above it follows that 
trace k = trace P and det k = det P. When these definitions are extended into each chart 
on A^ equation (fTTl ) shows that trace k G C'^{N) and det k G C°°{N). Moreover, if k is 
written as in equation ^, then 



tracer = —kY 



1 



At a point p £ N we say that k is invertible if (det k)\p ^ 0. If Id is the identity tensor 
Id G ^\{N), then writing Id as in equation © gives Id^^, = SlSi - Sldt. For / G C^{N) 
it follows that trace / Id = 6/. 



2.4. Decomposition of electromagnetic medium. At each point of a 4-manifold A^, an 
element of $7^2 (^) depends on 36 parameters. Pointwise, such (2) -tensors canonically de- 
compose into three linear subspaces. The motivation for this decomposition is that different 
components in the decomposition enter in different parts of electromagnetics. See IIHO031 
Section D. 1.3]. 
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Proposition 2.1. Let N be a A-manifold, and let 

Z = {KeQ^2{N) -.uAKiv) = k{u) Avforall u,v en^{N), 

tracer = 0}, 
W = {Ken'^2{N) :uAK{v) = -K{u)Avforallu,vGn'^{N)}, 
U = {fld£h\{N):f£C'^{N)}. 

Then 

(18) n\{N) = z ® w ® u, 

and pointwise, diniZ = 20, dimTy = 15 anddiuiU = 1. 

If we write a k E H^iN) as k =(i) k + C^k + (^k with (^k G Z, ^^V G M^, ^^k E C/, 
then we say that ^^k is the principal part, (^^k is the skewon part, ^^W is the axion part of 
K IIHO03L For a proof of Proposition 12. II as stated above, see HDahllaL and for further 
discussions, see IIRub02llHO03llFavl2L 

In 0^2 (-^) there is a canonical isomorphism 0^2 (-^) — ^ ^^2(-^) known as the Poincare 
isomorphism IIGre78l IFavl2ll . Let us first give a local definition. If k € ^'^2{N) on a 
4-manifold A^, we define k as the element k € J]^2(-^) defined as 

(19) K^g = —e^gab^cd^ 

when «; and k are written as in equation dlji. Equations (|7J-([8j imply that this assignment 
defines an element k G 0^2 (-^)- For k G ri^2(-^) we define k in the same way and we also 
have a canonical isomorphism $7^2 (-^) — ^ il^2(-^)- 

The next proposition collects results for k. In particular, part [[I}] states that ^ can be in- 
terpreted as a formal adjoint of k with respect to the wedge product for 2-forms. In con- 
sequence, the Poincare isomorphism is closely related to the decomposition in Proposition 
12.11 For example, k G ^1 (^) ^as only a principal part if and only if k = k and trace k = 0. 
For a further discussion, see nFavl21 . 

Proposition 2.2. Suppose N is a 4-manifold and k G 17^2 (-^)- 

(/) K is the unique k G 1^2 (A^) such that 

(20) k{u)Av = uAk{v) forall u,v en'^{N). 



(ii) /Id = fldforallj G C°°(iV). 
(///) K = K and ifr]£ r2^2(-^). ^'^^'^ ko rj = rj ok. 
(iv) trace k = trace k. 

(v) Ifu A k{u) = holds for all u G VL^{N) then k + k = 0. 



Proof. Part [(1)] follows by writing out both sides in equation (1201 ) in coordinates. Parts \(ii) 
and lf/Z/Jl foUow by part [[I}] Part lf/vJi is a direct computation. For part [fv}] we have 



_ 1 

u A {k + k){v) = - {{u + v) A k{u + v) — {u — v) a k{u — v)) 

for all n, u G ^"^{N), and the claim follows since the right hand side vanishes. D 
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If p is a twisted scalar tensor density of weight 1 on a 4-manifoId N and A, B £ Q,2iN) 
then we define pAiS>Ba.s the twisted tensor in Q'^2{N) defined as follows. If locally 

A = -2^'^^ A ^ and 5 = lB^^£, A ^ then 

(21) {pACi,B)i{ = persabA'^'B'^ 

when p A (g) i? is written as in equation ^. That pA^B transforms as an element in 
Q,'^2{N) follows by equation ([T]). Similarly when p is an untwisted scalar density we define 
pA^BG Q'^2{N) by equation (|2TI ). For both twisted and untwisted p we have identities 



(22) pA(g)B = pB®A, 

(23) {p'A®B)oK = pA(^{Bk), 



(24) Ko{pA(g)B) = p{Ak)(S)B, 

(25) {p^ B) o {p^ ^ A) = tiace{pB ^B) (pA^A). 

In Section l42l and in the proof of Theorem |5T| we will need the following lemma. 

Lemma 2.3. Suppose N is a A-manifold and k G Q'^2{N) is defined as 

(26) K = p {A(g)B + B0A)+ fid, 

where p is a scalar tensor density of weight 1, A,B G r22(-^) cind f € C°^{N). Then 
k\p = at a point p £ N implies that f\p = and p\p = or A\p = or B\p = 0. 

If K is written as in equation dUl and A, B are written as above, then equation (|26] | states 
that 






persab(A'^'B'^+A^^B''')+fld%. 



Proof. By restricting the analysis to p and introducing notation A^ = A^^^"^ and B^ = 
B^^^^, we obtain 

(27) 2p{A^B'^ + A'^B^) + fe^-^ = for all I,J £0. 

Setting I = J and summing implies that J2i&o pA^B^ = 0. Multiplying each equation in 
(IT71) by A^ B'^ and e^'^ and summing /, J yields two scalar equations. Eliminating / from 
these equations gives 

2\ 



V \/eO / V/eO / \i,j£0 




0, 
and the claim follows. D 

2.5. The Fresnel surface. Let k G Q'^2{N) on a 4-manifold N . If k is locally given by 
equation Q in coordinates {a;*}, let 

1 



/-90\ <^'J«:' _ _^i.aia2 azi a4,j bib2b^k b-ibuhal 

^^°> ^0 ~ 4g fei62 fe3b4 &5b6 taia2a3a4- 

If {P} are overlapping coordinates, then equations (l5j, (|7J and ([8]l imply that components 
'^Q satisfy the transformation rule 



(29) %^^^ 



dx^ 



cabcd ^^* ^^"^ '^^'^ ^^^ 

dx°- dx'^ dx'^ dx'^ 
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Thus components S^q-' define a twisted (g) -tensor density % on A^ of weiglit 1. The 
Tamm-Rubilar tensor density IIHO03I IRub02ll is the symmetric part of % and we denote 
this twisted tensor density by ^. In coordinates, ^^J'^' = ^j'-' ' , where parenthesis indi- 
cate that indices ijkl are symmetrised with scaling 1/4!. If locally ^ = ^idx^ it follows 

that '^'^^^iiijikii = %^^^iiijikiu and we call '^'^^^ iiijik^i the Fresnel polynomial. The 
Fresnel surface at a point p G A^ is defined as 

(30) Fp{K) = {eeT;(iv):^'^'='e.e,66 = o}. 

By equation (l29l) . the definition of Fp{K) does not depend on local coordinates. Let F{k) = 
UpeAf ^p('^) ^^ ^^^ disjoint union of all Fresnel surfaces. 

The Fresnel surface F{k) is a fundamental object when studying wave propagation in 
Maxwell's equations. Essentially, equation '^^^^''^iCjCk^i = in equation (l30l ) is a ten- 
sorial analogue to the dispersion equation that describes wave propagation in the geometric 
optics limit. Thus F{k) constrains possible wave speed(s) as a function of direction. In 
general the Fresnel surface Fp{K) is a fourth order polynomial surface in T*{N), so it can 
have multiple sheets and singular points IIOH04L 

There are various ways to derive the Fresnel surface; by studying a propagating weak sin- 
gularity MQFROOl IRub021 IHO03 1. using a geometric optics [Iti09[ |DahlTa] |. or as the char- 
acteristic polynomial of the full Maxwell's equations [SWWIO]. The tensorial description 
of the Fresnel surface is due to Y. Obukhov, T. Fukui and G. Rubilar HOFROOI . 

3. Results for skewon-free medium 

In this section we collect a number of results for twisted skewon-free tensors that we will 
need in the proof of Theorem 15. 1 1 

3.1. The normal form theorem by SchuUer et al. The normal form theorem for skewon- 
free medium by F Schuller, C. Witte and M. Wohlfarth HSWWIO I shows that there exists 
23 simple matrices such that any skewon-free medium can pointwise be transformed into 
one of these normal forms by a coordinate transformation plus, possibly, a conjugation by 
a Hodge operator. Next we formulate a slightly simplified version of this result that is 
sufficiently general for the proof of Theorem 15.11 Let us note that the original theorem 
in I.SWW10.I is formulated for area metrics. However, under mild assumptions these are 
essentially in one-to-one correspondence with skewon-free tensors in 0^2 (-^)- The below 
presentation is based on the reformulation in HDahllcl . 

Suppose L is an element in Q^{N) (g) ^i{N) on an n-manifold N. Then we can treat L 
as a pointwise linear map ^^{N) — )• Q^{N). By linear algebra, it follows that around 
each p £ N there are coordinates such that at p, components {Lj )ij is a matrix in Jordan 
normal form. Since there are only finitely many ways an n x n matrix can be decomposed 
into Jordan blocks, it follows that there are only a finite number of normal forms for L\p. 
It should be emphasised that the structure of the Jordan normal form is unstable under 
perturbations of the matrix. Hence, the normal form is in general only valid at one point. 
The normal form theorem in HSWWIOI is essentially an analogous result for skewon-free 
elements k in $7^2 (-^)- The difficulty in proving such a result is easy to understand. The 
matrix that represents k at a point is a 6 x 6 matrix. By a linear transformation in M^, we 
can transform this into an Jordan normal form, but such a transformation, a priori has 36 
degrees of freedom. On the other hand, for a coordinate transformation on A^, the Jacobian 
only has 16 degrees of freedom. It is therefore not obvious that coordinate transformations 
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have enough degrees of freedom to transform k into a normal form. See equation (ITtT i. For 
a further discussion, see HSWWlOIIDahllc J. 



The below theorem summarises the normal form theorem in IIS WW 101 specialised to the 
setting that we need here. Let us make three comments. First, the below theorem is formu- 
lated for twisted k G 0^2 (-^) instead of for area metrics in ISWWIOI (which are ordinary 
tensors) or untwisted k G f]^2(-^) in [DahllcJ . Second, the theorem contains the techni- 
cal assumption that k is invertible and the Fresnel surface has no 2-dimensional subspace. 
This greatly simplifies the result since it implies that there are only 7 possible normal forms 
and one does not need any conjugations by Hodge operators. These assumptions will also 
appear in Theorem 15.11 For a further discussion of these assumptions, see end of Section 
|5l Third, the reason the normal form theorem is useful can be seen from Proposition 12.11 
Namely, in arbitrary coordinates, a skewon-free k G il^2(-^) depends on 21 parameters. 
However, from Theorem 13. II we see that each normal form depends only on 2, 4 or 6 pa- 
rameters. This reduction of parameters will make the computer algebra feasible in Theorem 
EI 



The division into metaclasses in HSWWIOII is based on the Jordan block structure of the 
matrix representation of k at a point. Since this structure is unstable under perturbations, it 
can be difficult to determine the metaclass both in the numerical case and the symbolic case 
IILZW97L 



Theorem 3.1. Suppose N is a 4-manifold and k G il^2(-^)- tfp G N and 

(a) K. has no skewon part at p, 

(b) K is invertible at p, 

(c) the Fresnel surface Fp{K) does not contain a two dimensional vector subspace. 



Then there exists coordinates {x*}^^q around p such that the 6x6 matrix {kj)ij that 
represents k\p in these coordinates is one of the below matrices: 



Metaclass I: 



(31) 



/ai -/3i 
02 








as 



/32 

/33 













-/32 



Q2 





\ 



-/33 



"3 / 



• Metaclass II: 



(32) 





1 



-/3i 

Ql 



1 









a2 



/32 













/3i 



-/3i Ql 





\ 



-/32 



"2 / 
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(33) 



(34) 



(35) 



(36) 



Metaclass HI: 














/ 


ai 


-/3i 











\ 




/3i 


ai 
















1 





Ol 








-/3i 













Ol 


/3i 


1 










1 


-/3i 


Ol 





\ 





1 


/3i 








ai) 


Metaclass W: 
















^ai 








-/3i 





0\ 







02 








-/32 













03 








04 




/3i 








Ol 













/32 








02 







vo 





04 








03/ 


Metaclass V: 
















' ai 


-A 











0\ 




/3i 


Ol 
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Metaclass VI: 
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Metaclass VII: 
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(37) 



In each matrix the parameters satisfy 01,02,... G K, /3i , /32 , 

sgn /32 = • • • . 



G M\{0}anJsgn/3i 



Proof. Let ([/, x*) be coordinates around p, and let P = {nj)ij be the 6 x 6-matrix that rep- 
resents K at p in these coordinates. By treating [/ as a manifold with coordinates {x*}f^O' 
equation @ defines a tensor k G Q'^2{U). Since k is invertible at p and -Fp(k) has no 2- 
dimensional subspace, the Jordan normal form of P can not have a Jordan block of dimen- 
sion 2, . . . , 6 that corresponds to a real eigenvalue of P. For area metrics this is established 
in Lemma 5.1 in BSWWIOL (Or, for a translation to elements in Q\{U), see the proof 
of Theorem 2.1 in IIDahllbl .) In the terminology of HSWWIOI and MDahllbll this implies 
that kL is of Metaclasses I, . . ., VII. Hence Theorem 3.2 in BDahllcll (the restatement of 
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the normal form theorem in HSWWIOI ) impUes that around p, manifold U has a coordinate 
chart (C/, X*) such that at p, we have 

(38) TPT-^ = R, 

where T = {^^)ij is as in equation (IT6l) and R is one of the 6x6 matrices in equations 
(|3TI ) — (|37] ) for some parameters qi, a2, . . . G M and /3i, /32, . . . > 0. Since {U, x*) is a chart 
in A^ it follows that ([/, x*) is also a chart in N. Multiplying equation (1381 ) by sgn det ( ^ 1 

and comparing with equation dTTl ) shows that sgn det ( ^ ) R is the matrix that represents 



K G il.\{N) in coordinates {x*}?^g. If sgn det ( ^ ) = 1 or if i? is in Metaclasses I, 

IV, VI, VII, the claim follows. On the other hand, if sgn det ( ||^ j = — 1 and R is in 

Metaclasses II, III, V, it remains to prove that we can change the signs of the 1-entries in 
the normal forms by an orientation preserving coordinate transformation. Let {x*}?^q be 
coordinates determined by x* = Jjx-' for a suitable 4x4 matrix J = {Jj)ij- For Metaclass 
III a suitable Jacobian is {Jj)ij = diag(l, —1, —1, 1), and for Metaclass II and V a suitable 
Jacobian is 



J 



n 



/I 





0\ 





-1 





1 








\0 





1/ 



3.2. Non-birefringent medium. By a. pseudo-Riemann metric on a manifold A^ we mean 
a symmetric (2) -tensor g that is non-degenerate. If N is not connected we also assume that 
g has constant signature. By a Lorentz metric we mean a pseudo-Riemann metric on a 4- 

manifold with signature ( — 1-++) or (H ). Let (j be the isomorphisms (j : T*N — > TN, 

so that if locally g = Qijdx'^ ® dx^ then jj(aj(ix*) = aig^^ -A-. Using the (1-isomorphism we 
extend g to covectors by setting g{^, rf) = g{£}.,'if) when ^, 77 G T*{N). 

For a Lorentz metric g the light cone at a point p G A^ is defined as 

Np{g) = {i(^T;{N):g{i,i) = Q], 

and analogously to the Fresnel surface we define N{g) = IJ„gjv ^pid)- 

If 5 is a pseudo-Riemann metric on a 4-manifold N, then the Hodge star operator of g is 
defined as the *g G il^2(-^) such that if locally g = g-ijdx^ dx^ , and *g is written as in 
equation ^, then 



(39) {*g)i{ = ^Ad^\9"'9''eabrs, 

where det g = det gij and g^^ is the ijth entry of (g^ )~^. Then *g has only a principal part. 
See for example, IIHO031lFavl2] . Moreover, if ^r is a Lorentz metric and k = *g,we have 

(40) F{k) = N{g). 

Equation (l40l ) is the motivation for defining N{g) as a subset of the cotangent bundle. 

Definition 3.2. Suppose A^ is a 4-manifold and k G i7^2(-^)- Then k is non-birefringent if 
there exists a Lorentz metric g on N such that equation (l40l ) holds. 
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Thus, in non-birefringent medium, the Fresnel surface Fp{K) has only a single sheet, and 
there is only one signal speed in each direction. In non-birefringent medium it follows that 
propagation speed can not depend on polarisation. On N = M^, a specific example of a non- 
birefringent medium is k = y/-*g, where g is the Lorentz metric g = diag(— ^, 1, 1, 1) 



on M^. Then constitutive equation (fT2l ) models standard isotropic medium on M^ with per- 
mittivity e > and fi > 0. The next theorem gives the complete characterisation of all 
non-birefringent media with only a only a principal part. 

Theorem 3.3. Suppose N is a A-manifold. If k G 0^2 i^) satisfies '^•'k = 0, then the 
following conditions are equivalent: 

(i) ^'^'K = and k is non-birefringent. 
(ii) K^ = -fid for some function f G C°°{N) with / > 0. 
(Hi) there exists a Lorentz metric g and a non-vanishing function f G C°°{N) such that 

(41) K = f*g. 



Implication[(£}]^p]]was conjectured in 1999 by Y. Obukhov and F. Hehl IIOH99ilOFR00]| . 



Under some additional technical assumptions the imphcation was already proven in HOFROOl . 
However, the general case was only established in [FBll] by A. Favaro and L. Bergamin 
by a case by case analysis using the normal form theorem in [SWWIOJ . For an alterna- 
tive proof using a Grobner basis, see HDahllal and for similar results, see IILH04I IIti051 
IRRSllI and Section 13.31 below. Implication \(iii)\ =^ Vijl is a direct computation. In the 



setting of electromagnetics, implication \(ii)\ =^ \(iii)\ seems to first to have been derived by 
M. Schonberg |Rub02, SchVlJ. For further derivations and discussions, see liHO03i.iRub02[ 
IOFR00[IOH99llJad79l . 

When a general k G [7^2 (-^) on a 4-manifold N satisfies k^ = —/Id for a function 
/ G C°°{N) one says that k satisfies the closure condition. For physical motivation, see 
IIHO031 Section D.3.1]. For a study of more general closure relations, and in particular, for 
an analysis when k might have a skewon part, see IIFavl21lLBF121 . and Section |43] below. 

3.3. Medium with a double light cone. Since the Fresnel surface is a 4th order surface, 
the Fresnel surface can decompose into two distinct Lorentz null cones. In such medium 
differently polarised waves can propagate with different wave speeds. This is, for example, 
the case in uniaxial crystals like calcite l|BW99[ Section 15.3]. This motivates the next 
definition. 

Definition 3.4. Suppose A^ is a 4-manifold and k G 0.\{N). If p G N v/e say that the 
Fresnel surface Fp{K) decomposes into a double light cone if there exists Lorentz metrics 
g+ and g^ defined in a neighbourhood of p such that 

(42) Fp{K) = Np{g+) U Np{g_) 

md Np{g+)^Np{g^). 

If g, h are Lorentz metrics, then Np{g) C Np{h) implies that at p we have g = Ch for some 
C G M \ {0}. See for example IITou65ll . Thus, if k decomposes into a double light cone, 
then K is not non-birefringent. 

Under some assumptions, the next theorem gives the complete pointwise description of all 
medium tensors with a double light cone. The theorem generalises the result in BDahllbl 
to twisted tensors. 
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Theorem 3.5. Suppose N is a 4-manifold and k G il^2(-^)- Furthermore, suppose that at 
some p £ N 

(a) K. has no skewon part at p, 

(b) K is invertible at p, 

(c) the Fresnel surface Fp{K) factorises into a double light cone at p. 

Then exactly one of the below three possibilities holds: 

(i) Metaclass I. There are coordinates {x*}?^q around p such that the matrix (k'J)ij 
that represents k\p in these coordinates is given by equation (1311 ) for some ai , Q2 , 03 G 

M ancf /3i , /32 , /33 G M \ {0} with 

"2 = as, /32 = /^s, sgn /?! = sgn /32 = sgn /33 

and either ai 7^ 02 or /3i 7^ (32 or both inequalities hold, 
(ii) Metaclass 11. There are coordinates {x*}f^Q around p such that the matrix {kj)ij 
that represents k\p in these coordinates is given by equation (1321 ) for some ai,a2 G 
Rand/3i,f32 £R\{0} with 

ai = "2, /?! = /32- 
(Hi) Metaclass IV. There are coordinates {x*}f^Q around p such that the matrix (kj)ij 



that represents k\p in these coordinates is given by equation (1341 ) for some ai , Q2 , 03 , 04 G 

Rand/3i,f32 G M \ {0} w/f/i 

ai = "2, /3i = (32, 04, / 0, 03 / aj. 

Conversely, if k is defined by one of the above three possibilities, then the Fresnel surface 
of K decomposes into a double light cone at p. 

Proof. For k G Q,'^2{N) the result is proven in P Dahllbl Theorem 2.1] (up to a permutation 
of coordinates in Metaclass I). The generalisation to k G 17^2 (-^) follows by the same 
argument used to prove Theorem 13.11 The converse direction can be verified by computer 
algebra using the expUcit Lorentz metrics given in HDahllbl . D 



In Theorem 13.51 uniaxial medium is given by Metaclass I when ai = 02 = 03 = 0. The 
main conclusion of the theorem is that there are two (and only two) additional classes of 
medium where the Fresnel surface decomposes (Metaclasses II and IV). In all three classes, 
there are explicit formulas for the Lorentz metrics that factorise the Fresnel surface. For a 
further discussion of these metrics, see [Dahllb l. 

In Theorem 15.11 we will show that under suitable assumptions every skewon-free medium 
with a double light cone can be written as in equation (l43T l. This medium class is a special 
class of generalised Q-medium introduced by I. Lindell and H. Wallen in IILW02L For 
further discussions of this medium class, see IILW04I lFavl2l ILBF12I1 . 

Proposition 3.6. Suppose N is a A-manifold, g is a Lorentz metric, p is a twisted scalar 
density of weight 1, j4 G il2(-^) cind Ci G M \ {0} and C2 G M. Moreover, suppose 
K G il^2(-^) is defined as 

(43) K = Ci*g+pA(S)A + C2ld. 

Then k is skewon-free the following claims hold pointwise in N: 

(i) K is non-birefringent if and only if A = or p = 0. 
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(ii) K has a double light cone if and only if p ^ 0, A ^ and 
(44) det K / {Cf + Cif (C2 + ^ trace{pA (^ AU . 

Proof. We restrict the analysis to a point p £ N, and let {x*}?^g be coordinates around 
p such that the Lorentz metric has components g = ibdiag(— 1, 1, 1, 1) at p. For claim 



(i)\ let us note that the axion component of k does not influence the Fresnel polynomial. 
See for example IIHQ03I . Thus k is non-birefringent when ^4 = or p = 0. For the 
converse direction, suppose k is non-birefringent. Then Theorem 13.31 implies that {k — 
I trace K Id)^ = —Aid for some A > 0. Writing out the last equation and solving the 
associated Grobner basis equations (see HCLOOTIlDahl lal ) shows that A = or p = 0. For 
claim If //jl let us write A = ^A'^ -^ A A. Then the Fresnel polynomial at p is given by 

(45) ^'''''Ci^jUi = -c! [g'^Uj] {H'^^i^j) , 

where g'^ = {g~^)ij and H^^ = Cig'^ - IpA'^gabA''^ (see IILW02[ILBFT2]| '). Moreover, 

(46) det K = {Cf + C|) ^ {Cf + C^ + E + C2 trace(p A^A)), 

where i? G M is an expression that depends on p, Ci and A. We will not need the explicit 
expression for E. However, by computer algebra we see that the same E also appears in 
det H for matrix H = {W^)ij. Then equation (l46l ) yields 

2 



detH = -(cl + E-j{iTa.ce{pA®A)f\ 






(47) = - ( ,^2 , ^2N2 - ( ^2 + 7^ trace(pA 




If K has a double Ught cone, claim [fl}] implies that A ^ Q and p / 0. Moreover, by Propo- 
sition 1.5 in IDahllb] and since polynomials have a unique factorisation into irreducible 
factors [CLO07i Theorem 5 in Section 3.5], we have det H < and equation (l47l) implies 
inequality (l44l) for det k. Conversely, if the inequalities in claimlfZ/Jlare satisfied, then equa- 



tion (|47] ) shows that det H < 0, so g and H both have Lorentz signature at p. To complete 
the proof we need to show that there is no constant C G K\{0} such that g^^ = CWL Since 
A y^ and p 7^ 0, this follows by inspecting equations g" = CW^ for i = 0, . . . , 3. D 

4. Decomposable media 

In this section we first describe the class of decomposable medium introduced in fLBFlll. 
In particular, in Theorem 14.3 1 we describe the sufficient conditions derived in [LBF12J that 
imply that a medium is decomposable. In Theorem 15.11 these conditions will play a key 
role. In Section 1431 we will describe some results that suggest that condition [(7)] in Theorem 
I4.3l is a general factorisability condition for the Fresnel polynomial. Following IILBF12I1 we 
restrict the analysis to M^ so that we can work with plane waves. 

4.1. Plane waves in M^. We say that a tensor T on R^ is constant if there are global coor- 
dinates for M^ where components for T are constant. If we assume that many tensors are 
constant, we assume that they are constant with respect to the same choice of coordinates. 
Below we also use notation Q^{N, C) to denote the space of A;-forms on a manifold N with 
possibly complex coefficients. 
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Suppose K £ Q,\{R'^) is constant and F,G £ Cl'^{R'^) are defined as 

(48) F = Re{e**X}, G = Re{e"^Y}, 

where <I> is a function $ : M'^ — ;> M such that d^ is constant and non-zero, X,Y £ f!,'^{R^, C) 
are constant and not both zero. If F and G solve the sourceless Maxwell's equations we say 
that F and G is a plane wave. 

Proposition 4.1. Suppose k £ ri^2(I^^) is constant and $ is a function $ : M^ — t- M such 
that d^ is constant and non-zero. Moreover, suppose X, Y are constant 2-forms X,Y £ 
ri^(R^,C). If F and G are defined by equations (|48] ), then the following conditions are 
equivalent: 

(i) F and G is a plane wave. 

(ii) d^ £ F(k) and there exists a constant a £ ri^(R^,C) such that d^ A a 7^ 0, 
d^ A K{d^ A a) = and 

(49) X = d^Aa, 

(50) Y = K{d^Aa). 

Proof. Let ^ = d^.lf F and G is a plane wave then ^ / implies that 

(51) C^X = 0, ^AY = 0, Y = k{X). 

The first equation in equation (ISTT i implies that there exists a constant 1-form a £ Q^ (W^, C) 
such that X = ^ A a. It is clear that a and ^ A a are both non-zero, since otherwise 
X = Y = 0. Combining the latter two equations in equation (|5TI ) implies that 

(52) CAK(^Aa) = 0. 

Since this linear equation for a has a non-zero solution, it follows that ^ £ F{k). See 
for example, HOFROOI lRub02l IHQ031[Dahl lall . This completes the proof of implication IflJ 



=^ \(ii)\ For the converse implication it suffices to verify that equations (|48])-(|501) define a 
solution to Maxwell's equations. D 

4.2. Decomposable medium. The next definition and theorem are from IILBF121 . It is not 
known if the converse of Theorem 14. 3 l is also true [iLBF12il . 

Definition 4.2. Suppose k £ (7^2(1^^) is constant. Then we say that k is decomposable if 
there exist non-zero and constant A, B £ Q,2{^'^) such that if F, G is a plane wave solution 
to Maxwell's equations, then 

(53) F{A) = or F{B) = 0. 

Tiieorem 4.3. Suppose n £ il^2(I^^) is constant. Furthermore, suppose 



(i) there exists constant tensors A,B£ Q,2i^ ) cind a constant scalar density p of 
weight 1 such that 

(54) a Id +/3 (k + k) + 7K o K = p(A®B+^®A) 

for constants a, /3, 7 £ M and (3, 7 are not both zero, 
(ii) the right hand side in equation (1541 ) is non-zero. 



Then k is decomposable (and condition (1531) holds for the same A and B as in condition 
dH). 



Before the proof, let us note that by Lemma 12.31 the right hand side in equation (l54l) is 
non-zero if and only if ^, i? and p are all non-zero. 
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Proof. (Following IILBF12 II.) Suppose condition (l54l ) holds for some a, /3, 7, p, A, B. More- 
over, suppose F, G is an arbitrary plane wave for k as in equation (|48] ). To prove the claim 
we need to show that condition (l53l ) holds. Proposition 14. 1 l implies that Y = k{X) and 

XAX = 0, XAY = 0, YAX = 0, YAY = 0, 

whence equation (l20l ) implies that 

(55) = X A{ald+f3{K + K) + jKOK){X). 

Let {x*}?^Q be coordinates for R^ where all the aforementioned tensors are constant. Then 

= X Ap{A<S)B + B(^A){X) 

= X{A)X{B) pdx^ A dx^ A dx^ A dx^ . 

Here, the first equality follows by condition (l54l ) and (|55] ). and the latter equality follows 
by a computation in coordinates. Since A and B are real, it follows that F{A) = or 

F{B) = 0. n 

In Theorem l5.1l we will see that all the medium tensors in Theorem 13.51 are decomposable. 
In particular, uniaxial medium is decomposable. The next proposition shows that isotropic 
medium determined by a Hodge star operator is never decomposable. 

Proposition 4.4. Suppose k G i7^2(IK^) is defined as 

K = Ci*g + C2 Id, 

where Ci G M \ {0}, C2 G M and g is a constant indefinite pseudo-Riemann metric on M . 
Then k is not decomposable. 

Proof. Let us first assume that 5 is a Lorentz metric and let {x*}?^q be coordinates such 
that g = k diag(-l, 1, 1, 1) for some k e {-1, 1}. At G M^, it follows that 

fo{k) = {eGTo*(M"):-eo' + e? + e2 + ei = o}. 

For a contradiction, suppose k is decomposable. By Proposition 14. 1 1 there exists a non-zero 
and constant A, B £ ^20^'^) such that 

(56) {CAa){A) {CAa){B) = 
for all ^, a G T^{R^) that satisfy ^ G Fo{k) and 

(57) ^ A a / 0, ^ A k(^ A a) = 0. 

Let G is the subset G C Fo(k) \ {0} for which each coordinate belongs to {0, 1, \/2, \/3}. 
That is, one can think of G as a discretisation of Fq{k) in one quadrant of Tq*(M^). In 
total there are 19 such points, and for each ^ G G, we can find two linearly independent 
a G rg*(M^) such conditions ^} holds, cf BDahllall . Insisting that equation ^^ holds 
for all such ^ and a gives 19 x 2 = 38 second order polynomial equations for variables in 
A and B. Computing a Grobner basis for these equations and solving implies that either 
^ = or i? = 0. See JCLOOTI . Hence k is not decomposable. When g has signature 
( 1-+) the claim follows by repeating the above argument. D 
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4.3. Factorisability of the Fresnel polynomial. In what follows condition |[7}] in Theorem 
will play a key role. Let us therefore introduce the following definition. 



Definition 4.5. If k G il^2(I^^) is constant and satisfies condition [fl}] in Theorem 14.31 then 
we say that k is algebraically decomposable. 

In IILBF12L I. Lindell, L. Bergamin and A. Favaro showed that if k is algebraically de- 
composable (plus some additional assumptions), then the Fresnel polynomial of k always 
factorises into the product of two quadratic forms. In this section we summarise this result 
in Theorem 14.61 Moreover, we will see that for algebraically decomposable medium, the 
Fresnel polynomial seems to factorise even when the additional assumptions in Theorem 
14.61 are not satisfied. These results suggest (but do not prove) that the definition of alge- 
braically decomposable medium might be a sufficient condition for the Fresnel polynomial 
to factorise. 

Let us first note that the class of algebraically decomposable media contains a number 
medium classes as special cases. If k is purely skewon, then k + k = and k is al- 
gebraically decomposable. Also, if n satisfies the mixed closure condition k o k = A Id 
IILBF12irFavI2L then k is algebraically decomposable. If k has no skewon part, then k = k 
and the definition of algebraically decomposable medium simplifies. Thus, if k has no ske- 
won part and if k is a self-dual medium (so that a Id +/?«; + 7^^ = 0) MLinOSII . then k is 
algebraically decomposable. In particular, skewon-free medium that satisfies the closure 
condition k^ = A Id [HQ03J is algebraically decomposable. 

Equation (l54l ) that defines algebraically decomposable medium is a nonlinear equation in 
K. Suppose {x*}f^Q are coordinates for R^, P G M^^^ is the matrix P = {kJ)ij that 
represents k and A,B £ M.^ are the column vectors A = {A^)i and B = {B^)j that 
represent bivectors A and B with components as in Section 1241 Then equation (l54l ) reads 

(58) aE + I3{P*E + EP) + -yP^EP = 2p{AB^ + BA^), 

where A^ is the matrix transpose and E G M^^^ is the matrix E = {e^-^)ij. Numerically, 

E = { J ] , where and / are the zero and identity 3x3 matrices. When 7 7^ 0, equation 

(l58l) is structurally similar to an algebraic Riccati equation IIGLR05L 

The next theorem summarises the factorisation result from IILBF12I1 . but restated in the 
present setting. 

Theorem 4.6. If k £ 17^2 (I^^) is algebraically decomposable and a, /?, 7, p, A, B in equa- 
tion (1541) satisfy one of the below conditions: 

(i) 7 = 0, 
(//) 7 7^ 0, /3^ — 07 7^ and there exists a D £ ^^2(1^^) such that 

(59) D{-fK + pld) = ^tTace{pD(S>D)A + -fB. 

Then the Fresnel polynomial of k factorises into the product of two quadratic forms. 



Let us note that equation (I59b is a non-linear equation for D. A priori, the equation has real 
solutions, complex solutions, or no solutions for D. For a discussion of the last possibility, 
see below. Pointwise trace(/9 Z? (g) D) = holds if and only if D A D = or p = 0. 

Let us outline the argument in IILBF12II used to prove Theorem 14.61 Suppose ^\{M.^) is 



algebraically decomposable. If assumption (/) holds, then by rescaling we may assume that 
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/? = 1. Then, since k + k = 2{'-^'>k +(^) k), it follows that 

(60) aId+2(K-cj) = p{A(g)B + B<S)A) 

for some a € il^2(IR^) with only a skewon part. This gives an explicit representation 
formula for all k that satisfy condition (l54l ) with 7 = 0. Computing the Fresnel polynomial 
for K shows that it factorises into two quadratic forms. On the other hand, when assumption 
(ii)\ holds, then Theorem 14.71 in the below shows that equation (l54l ) transforms into r} o 



T] = Xld for some A / by a transformation similar to completing the square. Thus, 
to understand the structure of algebraically decomposable medium that satisfy assumption 
(ii)\ we only need to understand the simpler equation rj o r] = Xld with A / 0. In [LBF12II 



the latter equation is solved (see also IIFavl2ll ) using two explicit representation formulas 
similar to equation (|60] |. Using these representation formulas, the Fresnel polynomial can 
again be computed, and in both cases it factorises into a product of quadratic forms. 

The next theorem from fLBF12ll describes the transformation property of equation (|54] ) 
used in the proof of Theorem 14.61 The proof is a direct computation using identities (l22l) - 
(I25] ). For a general discussion of transformation properties for the matrix algebraic Riccati 
equation, see IICPL10[ILRT21 . 



Theorem 4.7. Suppose k £ 1)^2 (1^^) is algebraically decomposable such that equation 
(1541) holds with 'y ^ 0. If, moreover, there exists a D £ $72(1^^) such that equation (1591 ) 
holds, then rj G 0^2(1^'') defined as 



(61) rj = -iK- pD® A + pid 
satisfies 

(62) 7? o 7? = (/32 - 07) Id . 

Suppose K is algebraically decomposable such that equation (l54l ) holds with 7 7^ and 
/3^ — 07 = 0. Now we can not use Theorem 14.61 do decise whether the Fresnel polyno- 
mial factorises. However, by computer algebra we can find explicit examples of medium 
tensors with the above properties. Preliminary computer algebra experiments using such 
expressions suggest that the Fresnel polynomial always seems to factorise when the above 
assumptions are met. However, the factorisation seems be qualitatively different. Condition 
/3^ — 07 = seems to imply a linear factor in the Fresnel polynomial. For example, the 
Fresnel polynomial can factorise into the product of irreducible 1st and 3rd order polyno- 
mials. On the other hand, suppose k is algebraically decomposable such that equation (l54l ) 
holds with 7 / 0, /3^ — 07 7^ and equation ( [59l ) has no real solution for D. Now we can 
neither use Theorem l4.6l do decise whether the Fresnel polynomial factorises, but we may 
again construct explicit examples of medium tensors with the above properties. Using these 
expressions, preliminary computer algebra experiments suggest that the Fresnel polynomial 
also seems to factorise in this case. In conclusion, these initial observations together with 
Theorem 14.61 suggest that the definition of algebraically decomposable medium could be a 
sufficient condition for the Fresnel polynomial to factorise. 

Lastly, let us note that algebraic Riccati equations, and more generally, quadratic matrix 
equations, appear in a number of fields. In view of Theorem 14.61 and equation (|58] ). it is, 
however, interesting to note that quadratic matrix equations appear in the study of polyno- 
mial factorisation in one variable [BG05 1. Differential Riccati equations also appear in the 
problem of factoring linear partial differential operators of second and third order [GS041 . 
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5. Characterisation and representation of media 
with a double light cone 

Theorem 5.1. Suppose N is a A-manifold, and k G VL^2{^) is skewon-free and invertible 
at a point p & N. Then the following conditions are equivalent: 

(i) The Fresnel surface of k decomposes into a double light cone at p. 
(ii) K satisfies conditions: 

(a) the Fresnel surface Fp{n) C T*{N) does not contain a two-dimensional vec- 
tor subspace. 

(b) there are A,Bg O.'^(N) and a tensor density p of weight 1 such that at p we 
have 

(63) {K + iildf = -\ld+p(A®B + BCdA) 

for some p G C'°^(iV) and A G C°°{N). Moreover, A,B, p y^ and X > at 
p. 
(Hi) Around p there is a locally defined Lorentz metric g, a locally defined non-zero 
twisted scalar density p of weight 1, an A £ 0,2i^) that is non-zero at p, and 
constants Ci G M \ {0} and C2 G M such that at p, 

(64) K = Ci*g+pA(g)A + C2ld, 
and K satisfies inequality (1441) at p. 

As described in the introduction, the above theorem is the main result of this paper. A 
discussion of the theorem is postponed to the end of this section. 

In the Theorem 15. II we will use the computer algebra technique of Grobner bases IICLQ07II 
to eliminate variables from polynomial equations. This technique was also used in MDahllbll . 

Let C[ui, . . . ,UAr] the ring of complex coefficient polynomials C^ — )• C in variables 
ui, . . . ,un- For polynomials ri, . . . ,rfc G C[ni, . . . ,un], let 

k 
(ri,...,rfc) = (^fin : ft G C[ui, . . . ,un]} 

i=l 

be the the ideal generated by ri, . . . ,rk. Suppose V C C^ is the solution set to polynomial 
equations pi = 0, . . . ,pM = where pi G C[ui, . . . , un]. If / is the ideal generated by 
Pi, ■ ■ ■ ,Pm, the elimination ideals are the ideals defined as 

h = InC[uk+i,...,UN], k£{0,...,N-l}. 

Thus, if (ui, . . . ,un) G V then by IICLO07[ Proposition 9, Section 2.5] it follows that 
p{uk+i, ■ ■ ■ ,un) = for any p G Ik, and I^ contain polynomial consequences of the 
original equations that only depend on variables u^^i, . . . , U]\j. Using Grobner basis, one 
can explicitly compute I^ IICLQ07i Theorem 2 in Section 3.1]. In the below proof this has 
been done with the built-in Mathematica routine 'GroebnerBasis'. The same technique of 
eliminating variables was also a key part of the proof of Theorem [33] in BDahllbi 



Proof. Let us first prove implication \(i)\ =^ \(ii)\ By MDahllb l Proposition 1.3] condition 



[(7)] implies that Fp{K) has no two dimensional subspace. By Theorem 13.51 we only need to 
check three medium classes. 

Metaclass I. If k|p is in Metaclass I, then k can be written as in equation (|3TI ) with con- 
ditions on the parameters given by Theorem 13.51 Suppose ai = 02- Then Theorem 
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implies that /3i / /32- Let p = i(/3| - /3f), /x = -ai, A = /3|. Moreover, let A and B be 



bivectors defined as A 



,A'^ 



dx' 



A 



dx^ 



and similarly for B, with coefficients 



(65) 



(^A^r 



V 



/O 1 o\ 




/O o\ 






, iB%j = 



1 


V V 




V 0/ 



where subdiagonal terms are determined by antisymmetry. For these parameters, computer 
algebra shows that equation (1631 ) holds. On the other hand, if ai 7^ 02. suitable parameters 
are 



P 



8{ai - a2)h ' 



M 



-"2, 



and 



^A^r 



fO 2(ai - 02)/?! 




u 



A = /32^ 








V 



(ai - 02)^ - /3f + /5i + ^/^ 




where 



a = ((ai - a2)^ + (/3i - (32?) ((ai - 02)=^ + (/3i + (32?) . 
Bivector B is defined by the same formula as for A, but by replacing ^/a with —^. 



Metaclass II. If k\p is in Metaclass II, then k can be written as in equation (1321 ) with condi- 
tions on the parameters given by Theorem 13. 5 1 Suitable parameters are p = /3i/2, ^u = —a\, 
A = /3f and 



(66) 





/O 1 1 o\ 


[A'')^, = 








V 0/ 



{B^^) 



u 



/O 1 -1 0^ 





V 0/ 



Metaclass IV. If k|p is of Metaclass IV, then k can be written as in equation (l34l ) with 
conditions on the parameters given by Theorem 13.51 If a\ 7^ 03, then suitable parameters 
are 



1 



8(03 — ai)a4 



/^ 



-Qi, A = /3i 



and 



(A^^-) 



/O (ai - 03)^ + al + /3? + ^\ 

2(03 — ai)a4 



V ) 



where 



a=[a\- (03 - oixf) + /?? (204 + f\ + 2(ai - 02)^) 
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and B is defined as in Metaclass I. On the other hand, if ai 

\{l3l + al), /i = -as, A = /3f and 



as, then suitable parameters 



are p 



(67) 



{A^^) 



«j 



/O o\ 




/O l\ 


1 



, iB%j = 






V V 




^ 0/ 



This completes the proof of implication (i) =^ (ii) 



For the converse implication \(ii)\ =^ (i) suppose that k satisfies the conditions in \(ii)\ By 



Theorem 13. II we may assume that there are coordinates {x*}?^q around p such that at p, 
tensor k is given by one of the matrices in equations (|3TI)-(|37]) for some parameters as in 
Theorem 13. II Let us consider each of the seven cases separately. 

Metaclass I. If k\p is in Metaclass I, then there are coordinates {x*}?^q around p such that 
K is given by equation (|3T]) . By scaling A and B we may assume that p\p = 1. Moreover, 
writing out equation (1631 ) and eliminating variables in A and B using a Grobner basis (see 
above) yields equations that only involve A, fi and the parameters in k. The rest of the 
argument is divided into three subcases: 

Case 1. If /3i = /32 = /^a the Grobner basis equations imply that A = /3^ and 



(68) 
(69) 

(70) 



(a2 + ^)(a3 + )u) = 0, 
{ai + p){a3 + fi) = 0, 
(ai + /i)(a2 + /i) = 0. 



It follows that ai , a2 , as can not be all distinct, and by a coordinate change, we may assume 
that a2 = as- If ai = a2 = as, equation (l68l ) implies that /i = — ai. Then equation (|3T]) 
implies that k = — /3i *g +ai Id at p, where g is the Hodge star operator for the locally 
defined Lorentz metric g = diag(— 1, 1,1,1). Then equation (l63l) implies that p{A® B + 
B ® A) = 0. Since this contradicts Lemma [231 we have ai 7^ a2 and k has a double light 
cone at p by Theorem [33] 

Case 2. If exactly two of /3i , /32 , /3s coincide, then after a coordinate change we may assume 
that /3i 7^ /32 = /3s. Then the Grobner basis equations imply that either A = /3^ or A = /3|. 
If A = Pf, the Grobner basis equations imply that ai = 02 = as and /3i = /32 = /33. 
We may therefore assume that A = /3|. Then the Grobner basis equations imply that p = 
—02 = —as, and k has a double light cone at p by Theorem 13. 5 1 

Case 3. If all Pi, 132, Ps are all distinct, then the Grobner basis equations imply that 



(/3|-A)(/3|-A)(ai+^) 
(/3f-A)(/3|-A)(a2 + ^) 
(/3f-A)(/3i-A)(a3 + /i) 



(/3f-A)(/3| 



A)(/3| 



A) 



0, 
0, 
0, 
0. 



These equations imply that we must have A = /3f and p = —ai for some i € {1, 2, 3}. If 
1 = 1 the Grobner basis equations imply that ai = a2 = as and /?i = (32- This contradicts 
the assumption that all /3j are distinct. Similarly, i = 2 and i = 3 lead to contradictions, and 
Case 3 is not possible. 

Metaclass II. If k\p is in Metaclass II, there are coordinates {a;*}?^Q around p such that k is 
given by equation (l32l ). Writing out equation (|63] ) and eliminating variables as in Metaclass 
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I gives equations that only involve variables A, fi and the variables in k. Solving these 
equations give 

Ij, = -a2, A = /3|, /3i=/32, ai = a2, 



and K has a double light cone at p by Theorem [ 

Metaclass III. If k\p is in Metaclass III, there are coordinates {x*}f^Q around p such that k 
is given by equation (|33] ). Eliminating variables as in Metaclass I implies that /3i = 0. Thus 
k\p can not be in Metaclass III. 

Metaclass IV. If k\p is in Metaclass IV, there are coordinates {x*}f^Q around p such that k 
is given by equation (l34l ). We have 04 7^ since otherwise span{dx^\p, dx'^\p} C Fp^n). 
Moreover, since k is invertible at p it follows that 0:3 7^ Q4. Writing out equation (l63T l. 
eUminating variables as in Metaclass I, and solving implies that 

A = /3i, Pi = 132, /i = -Qi, ai = a2, 
and K has a double light cone at p by Theorem 13.51 

Metaclass V. If k\p is in Metaclass V, there are coordinates {x*}^^q around p such that k is 
given by equation (|35T ). We may assume that 03 / 0, since otherwise span{(ix*|p}f^^ C 
Fp{K). Eliminating variables as in Metaclass I, and solving implies the contradiction A + 
a\ = 0. Since A > it follows that k\p can not be in Metaclass V. 

Metaclass VI. If K,\p is in Metaclass VI, there are coordinates {x*}?^q around p such that n 
is given by equation (l36l ). Eliminating variables as in Metaclass I implies that 

{\ + al + {a■i + ^lf){\ + {a2-ai + ^lf){\ + {a2 + ai + |If) = 0. 

Since A > 0, it follows that k\p can not be in Metaclass VI. 

Metaclass VII. If k\p is in Metaclass VII, there are coordinates {x*}?^q around p such that 
n is given by equation (|37] ). Eliminating variables as in Metaclass I and solving impUes that 

3 

fc=i 

Since A > 0, it follows that k\p can not be in Metaclass VII. This completes the proof of 



implication (ii) =^ (i) 



Implication \(iii)\ =^ [(?}]is a restatement of Proposition 13.61 To prove implication [fj)] =^ |fn/j 



we proceed as in implication [fj}] =^ \(ii}\ and by Theorem 13.51 we only need to check three 
medium classes. Also, by Proposition 13.61 we do not need to prove inequality (l44l) since it 
follows form the other conditions in |f /// j| when [(7}]holds. 

Metaclass I. If k\p is in Metaclass I, there are coordinates {x*}?^q around p such that k is 
given by equation (|3TI ) with conditions on the parameters given by Theorem 13.51 Suppose 

Qi = ao. Let Ci = . ^ , Co = ao, ^ = ^ and in coordinates jxH, let p be 

*-^|detg| PI 

defined by p = (/3| - j3l)/{2l3i). Then equation dMll holds when A = ^B'^£^A ^ when 
coefficients B^^ are as in equation (1651 ) and (7 is the Lorentz metric g = gijdx'' dx^ with 
coefficients 

( ( ^ ^ 

(71) (ffij)ii = (^'^S( 1,-1,-^,--^ 
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On the other hand, suppose ai 7^ 02- Let 'I' be one of the two roots to the quadratic equation 



(72) 



1 



■^^ -D^^ + h 



0, 



where 1)3 is defined as in MDahllbl Theorem 2. 1 f /j] 



D-, 



hh 



Since sgn/3i = sgn/32, the discriminant of equation (1721 ) is strictly positive. Thus ^ G 

R \ {0} and sgn ^ = sgn /3i. Let H G M be defined as 



Since ai / 02 we see that ^ = -g^- is not a solution to equation (1721 ) whence H / 0. Let 
Ci, C2 be as in the ai = 02 case and let p = sgn H. Then equation (l64l) holds when g' is 
the Lorentz metric given by equation (1711 and A = ^A^^ -§^ A A- is given by 

/O v^ 




{A'- 



)ij 









V 






ai~-Q2 





Metaclass 11. If k\p is in Metaclass H, there are coordinates {x^j^^q around p such that k is 
given by equation (l32l ) with conditions on the parameters given by Theorem [33] Let Ci = 

C2 = ai and p = 1/2. Then equation dMll holds when A = ^ A*-? ^ A ^ is 

and ^r is the Lorentz metric g = gijdx^ ® dx^ with coefficients 

-1 



I3ly/\detg\ 

as in equation 



(73) 



{9i. 



vnj 



/-I 








/3i\ 





-/3i 














-A 





V/3i 








0/ 



Metaclass IV. If k|p is in Metaclass IV, there are coordinates {x*}^^g around p such that k 



is given by equation 
Qi = Q3. Let Ci = 



with conditions on the parameters given by Theorem 13. 5 1 Suppose 
,(72 = 01,^ = a4,//3i and p = [al + /3f )/(2a4). Then 



*\/|detg| 



equation (1641 ) holds when yl = ^5*-^ A- A A- when i?*^ are as in equation (l671 i and (/ is the 
Lorentz metric (7 = gijdx'^ dx^ with coefficients 

(74) (5ij)ij = (diag(l,^,^,-l))-^ 

On the other hand, suppose ai 7^ 03. Let ^ be one of the two roots to the quadratic equation 

(75) ^2_^L>i^-l = 0, 
where (see BDahllbl Theorem 2.1 (Hi)]), 



D, 



(02 - 03)^ + Pi 



Ua 



hOiA 



Then 'I' € M \ {0} and since ai / ^3 equation (jTSj) imphes that ^ / |^. Thus H € M \ {0} 
when 

H = 1 (a^ _ /3ixi/) . 
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Let Ci, C2 be as in the ai = 03 case and let p = sgnS. Then equation (1641 ) holds when 

g is the Lorentz metric in equation (174] ! and A is the bivector A 

coefficients 



1 Aij -d_ 
2^ dx' 



A g|j with 



(^^r 



^J 




This completes the proof of implication (i) =^ (Hi) 



D 



Let us first emphasise that the conditions in Theorem 15.11 are written analogously to the 
conditions in Theorem 13.31 In each theorem, condition ^I)\ is the dynamical description 
of the medium, condition \(ii)\ is a characterisation of the medium and condition \(iii)\ is a 
general representation formula. Let us also emphasise that in suitable limits, condition 

—Aid in Theorem 13.31 and 



in Theorem 15.11 reduces to the closure condition k^ 
representation formula (l64l in Theorem 15.1 [ reduces to k 



f*g in Theorem 13.31 Let us 



also emphasise that in both theorems, all conditions are tensorial, and do not depend on 
coordinate expressions. A difference between the theorems is that Theorem 13. 3 1 is a global 
result, while Theorem 15. H is a pointwise result. 

All the mediums in Theorem 15.11 satisfy the technical assumptions in Theorem 14.61 with 
either D = Aor D = B when A and B are as in equation (1631 ). 

As described in the introduction, condition (ii) in Theorem 15. 1 l is a slight strengthening of 
the conditions derived in IILBFI21 (see Theorem 14.3 1 in the above). Representation formula 
(l64l) in Theorem 15. II is also adapted from [LBF12|. For constant medium tensors on M^, 
Theorem 15.11 implies that if k is invertible, skewon-free and has a double light cone, then 
K is algebraically decomposable, and hence decomposable by IILBFI2I (see Theorem l4.3l ). 
In this setting. Theorem 15.11 explicitly shows that the behaviour of signal-speed imposes 
a constraint on the behaviour of polarisation. This can be seen as somewhat unexpected. 
However, the explanation is that polarisation and signal speeds are not independent for a 
propagating wave, but constrained by equation (l52l ). For a further discussion, see lIDahllal . 
It is also instructive to note that condition (|63] ) is a second order polynomial constraint on 
the coefficients in k, but the definition of a double light cone involves the Fresnel surface, 
which is a constraint involving third order polynomials of the coefficients in k. The same 



phenomenon appears in equivalence (i) <^\(u)\m Theorem 



Part of condition \(ii)\ is condition (a), that states that the Fresnel surface of k contains 
no two dimensional subspace. Let us describe five results where this condition also ap- 
pears. First, if the Fresnel surface of a k G Q,'^2{N) can be written as Fp^n) = {^ G 
T*{N) : (g(C;C))^ = 0} for a pseudo-Riemann metric g, then condition (a) is satisfied if 

and only if g has signature ( 1-+). This follows by a result of J. Montaldi HMonOTI . 

For example, ii g = diag(— 1, —1, 1, 1), then Fp{K) contains the 2-dimensional subspace 
span{^|(y + ^p-, -A- + ■^}- Second, one can prove that condition (a) is always satisfied 
if K decomposes into a double light cone (Proposition 1.3 in BDahllbll ). Third, in matter 
dynamics systems, condition (a) can be motivated by the behaviour of energy [RRSlll. 
In the terminology of LRRSllJ . condition (a) can be replaced by the stronger condition 
that K is bihyperbolic. Fourth, condition (a) also appears in the study of the well posed- 
ness of Maxwell's equations as an initial value problem IS WWIOH . Lastly, in the normal 
form representation of skewon-free medium tensors in [iSWWlOH . condition (a) simplifies 
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the representation since the condition excludes all but the first 7 coordinate representations. 
See nSWWlO l and Section O in the above. 

When equivalence holds in Theorem 15.11 there does not seem to be a simple relation be- 
tween parameters Ci, C2, p, A, g in equation (l64l ) and parameters /i, A, />, A, B in equation 
(|63] ). However, if equation (l64l ) holds for an A such that ^ A vl = (that is, A is decom- 
posable or simple ICoh05 , p. 185]), then equation (|63] ) holds for parameters 



P 



-Co, A 



-cl 



B = A{*g 



Using a Grobner basis argument one can show that the tensor k defined by equation (|3TI ) 
when /3i = /32 = /33 = 1, ai = 1 and a2 = as = 2 is invertible and has a double light 
cone. However, it can not be written as in equation (|64] | for an A such that A f\ A = Q. 
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